Oscillations in the expression of a self-repressed gene 
induced by a slow transcriptional dynamics. 
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We revisit the dynamics of a gene repressed by its own protein in the case where the transcription 
rate does not adapt instantaneously to protein concentration but is a dynamical variable. We derive 
analytical criteria for the appearance of sustained oscillations and find that they require degradation 
mechanisms much less nonlinear than for infinitely fast regulation. Deterministic predictions are 
confirmed by stochastic simulations of this minimal genetic oscillator. 

PACS numbers: 87.18.-h 87.18. Vf 87.16.Yc 82.40.Bj 
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Networks of genes interacting via regulatory proteins 
modulating their activities are highly nonlinear systems 
which display a variety of dynamical behaviors [l|, 0, Q . 
Their modeling has generally assumed that gene activa- 
tion is fast compared to other processes so that transcrip- 
tion rate reacts instantaneously to protein concentration. 
However, transcription is a complex process In the 
last years, it has been increasingly recognized that gene 
activity fluctuations can be slow and that this can affect 
the behavior of gene regulatory networks. In particular, 
slow transcriptional bursting and transcriptional mem- 
ory have been observed experimentally [a H, 0]. The- 
oretically, it has been shown that slow activation dy- 
namics can lead to bursts in expression [§], induce 
bistability 10] or modify the flipping rate of a genetic 
switch 

In this Letter, we show that slow promoter dynamics 
can also lead to oscillations by investigating how tran- 
scriptional dynamics modifies the behavior of a single 




gene repressed by its own protein [l^, 0, on 0, m m, 

HI 13] • This old problem of theoretical biology has been 
recently revived by the study of the Hesl gene involved 
in the somite clock [2l|. The usual view is that oscilla- 
tions appear in this genetic circuit only when additional 
steps are inserted in the feedback loop 0, [3]. In the 
Goodwin and Bliss oscillators HI, UdUM, the gene pro- 
tein catalyzes synthesis of the actual repressor. In early 
circadian models, transport of the repressor into the nu- 
cleus [3, E3 is a key oscillatory ingredient. In fact, the 
mere introduction of a time delay in the one-gene circuit 
model (accounting for protein trans por t or more gener- 
ally a cascade of intermediate steps [22]) can destabilize 
it [H, [3, [13, H3]. As will be of particular interest here, 
oscillations may also be induced by strongly nonlinear 



degradation mechanisms [231 ] . 

For simplicity, we study the case of an elementary 
kinetic equation describing regulation through protein- 
DNA binding 24], which we however prefer to view as 
a minimal description of transcriptional memory in more 
complex mechanisms. We derive an analytical expression 
of the oscillation threshold, and show that when the gene 
response time is appropriately tuned, the one-gene circuit 
can be destabilized (and oscillations induced) by degra- 
dation mechanisms much less nonlinear than for infinitely 
fast regulation. This result provides new insights into the 
interplay of nonlinearity and time delay. Stochastic simu- 
lations confirm that the main results of our analysis carry 
over to low copy number situations. 

Our study is based on the following three-variable 
model describing the genetic circuit represented in Fig. [I] 



G = 6 (1 - G) - a C(P)G 
P = nG + /3 M-6 P F(P) 
M = mo + A G - 5 M H(M) 
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where G, P and M represent gene activity, protein and 
RNA copy numbers. Eq. (fTaf) formally describes the ki- 
netics of protein-DNA binding at rate cxq and unbind- 
ing at rate 9q [24[. Possible cooperativity effects are 
taken into account via the function C(P) and the num- 
ber n of proteins binding DNA. Single-protein regula- 
tion corresponds to C(P) = P, n = 1. More gener- 
ally, Eq ([lap is a minimal model for the dynamics of 
an effective gene activity G slowly relaxing towards an 
equilibrium value given by the gene regulation function 
G = 1/[1 + C(P)/C(P )}, with P the half-expression 
threshold. Such a model can be obtained as the leading 
approximation of a mechanistic model of transcription in- 
cluding all processes concurring to gene expression (e.g., 
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FIG. 1: Reaction diagram of the self-regulated gene circuit. 



protein-DNA interaction, formation of open complex, 
polymerase recruitment, chromatin remodeling,...) when 
there is a dominant limiting step. Eq. (jTaj) causes gene 
activity G to lag fast changes in protein level, and plays a 
dynamical role similar to explicit time delays [3, [2(| 
or to transport equations [12] . In Eq. (flbf . the three 
terms correspond to binding/unbinding, translation and 
degradation. Eq. (flc|) describes transcription at rate 
fio + XqG and RNA degradation. In order to understand 
how oscillations can be induced by tuning protein and 
RNA degradation, we derive the oscillation criterion for 
arbitrary degradation functions F(P) and H(M) with 
unit derivative at zero, 5p and 5m being the low-copy- 
number degradation rates. 

By renormalizing time, variables, parameters, cooper- 
ativity and degradation functions according to: 
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cW = ^), /w = ffi), ft(m) = 5M, (2d) 
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Eqs. ([T]) can be rewritten in dimensionless form 

g' = e[l-g(l + c(p))} (3a) 
p' = na [1 - 5 (1 + c(p))] + S[m - f(p)] (3b) 
m = p + Xg — h(m) (3c) 

where x' = dx/dt'. When / and g are monotonous 
and h[f(oo)] > p, model (J3J) has a single steady state 
(g*,p*,TO*) satisfying the fixed point equations: 



9* 



l + c(p.) 



j> / \ h(m*) - \i 
,m* = f(p*),g* = t- • (4) 



Note that the steady state depends only on parameters 
A and fj, as well as on functions c, / and h, whereas pa- 
rameters 9, a, 5 control time scales. The behavior of the 



degradation and cooperativity functions in the neighbor- 
hood of the steady state is described by the slopes 
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In the case of linear degradation [f(p) = p, h(m) = m], 
we have u = s = 1. Small or even negative values of 
the slopes s and u generally denote strongly nonlinear 
degradation mechanisms 0, |23| , including saturation. 

To assess whether Eqs. Q can display oscillations, we 
have searched for parameter values where the fixed point 
specified by ([4]) loses stability to a periodic solution via 
a Hopf bifurcation (i.e., a pair of conjugate eigenvalues 
of the linearized problem cross the imaginary axis). For 
simplicity, we assume perfect repression (jj, = 0) and a 
large threshold Po {ol ~ 0). Under this approximation, 
the Routh-Hiirwitz stability criterion |25( indicates that 
a Hopf bifurcation occurs when the quantity 



H = a + (SXvgl +<t 2 )t + jar 2 



(5) 



crosses zero to become negative, where r = g*/9 is the 
gene response time, and the sum a = 5s + u and prod- 
uct 7 = 5su are symmetric functions of degradation 
rates 5s and u. At bifurcation, the oscillation period 
is t osc = 2ivy/ t I '{a + yr) where a and 7 satisfy Ti. = 0. 
Cooperativity essentially changes feedback strength from 
5X to 5Xv, thus we assume for simplicity single-protein 
regulation [c(p) — p, v — 1] thereafter. 

Eq. ([5]) shows that a strong feedback destabilizes the 
system while high degradation rates (large a and 7) tend 
to stabilize it. In the single protein case, H > when 
both protein and RNA are linearly degraded and no 
oscillations occur. Conversely, when protein and RNA 
degradations are completely saturated (s = u = a = 0), 
H = —5Xg 2 r < 0, indicating that oscillations then ap- 
pear systematically. The behavior in intermediate cases 
depends on the value of the response time r. 

In the classical case r = 0, Ti = a and oscillations 
appear only for a < 0. It is indeed known that negative 
effective degradation rates can lead to oscillations [2, 23]. 
We thus restrict ourselves to showing that at finite r, os- 
cillations can occur for u, s > 0. More precisely, we want 
to understand how oscillations arise away from the sat- 
urated cases u = or s = 0. To this end, we use the 
geometric slope average v = ^/us = as an index 

(v = 1 in the linear case), seeking to determine the max- 
imum value of v at which oscillations can be observed, 
and for which values of r this extremum is achieved. 

The quantities a and 7 play complementary roles. The 
total degradation rate a controls instability onset for 
small to moderate r. Moreover, Eq. §5§ indicates that 
when a, 7 > 0, a necessary condition for oscillations is 



a < a c 



g*V5X 



>V2 



(6) 
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where t sw is the time during which a fully active gene 
synthesizes the amount of protein corresponding to half- 
repression threshold. The degradation rate product 7 
is relevant only for large r, blocking oscillations if it is 
too large. In particular, 7 = guarantees the onset of 
oscillations for sufficiently large r whenever © holds. 
For u, s > and a given value of a, 7 can take any 
value between a 2 /4 and depending on whether the two 
degradation rates are equal or completely unbalanced, 
one being equal to zero and the other to a. We use below 
e = 2^/7/cr S [0, 1] as a balance indicator. 

Remarkably, we note that under the rescaling 
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the oscillation condition can be rewritten without explicit 
parameter dependence: 



W e (E,T) = S x 



[e 2 E 2 , 
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<0 (8) 



and defines a series of curves £ e (T) such that given a 
balance index e and a response time T, oscillations are 
found for £ < £ e (T). Fig.|lshows the limit curves £ X (T) 
and So(T) which are important to understand the bifur- 
cation diagram: regardless of the value of e, the circuit 
always (resp., never) oscillates when £ < £j.(T) [resp., 
£ > £q(T)]. To support our analysis, we have searched 
the parameter space of Eqs. ([3]) for oscillatory behavior 
for a,fj,^0, assuming for definiteness allosteric protein 
degradation and Michaelis-Menten RNA degradation: 



px (a+p/n) 
a + 2a(p/n) + (p/n) 2 



h(m) 



X - 



(9) 



Points in the (£, T) plane associated with oscillating pa- 
rameter sets are shown as black dots in Fig. [5] Agreement 
is excellent: all dots are below the £o(P) curve and the 
few significantly above £i(T) have one small degradation 
rate. We are thus confident that our analysis allows us 
to understand the behavior of Eqs. (|3|). 

Two regions can be distinguished in Fig.O For T < 1, 
the instability threshold £ e (T) is practically indepen- 
dent of e and increases rapidly with T. For small T, 
£ e (T) ~ T (thus, the oscillation criterion is a < <j 2 t). 
In the T > 1 region, £ e (T) reaches its maximum value 
£m(e) = VVeTT at T = T m {e) = 2y/e+l/e, and then 
decreases as T^ 1 for T — > 00, except for e = where 
it monotonously increases towards £ = 1. At fixed £, 
oscillations are thus found in a finite range of T, which 
widens gradually, and is eventually infinite, as e — ► 0. 

In the oscillation region, the index v ~ e£ measur- 
ing distance from saturated degradation reaches its max- 

~ for T = T opt = 2y/2 and 



opt 



imal value 

e = 1, at the maximum of the £i(T') curve. Our analysis 
thus unveils a resonance phenomenon in the dynamics of 




FIG. 2: Bifurcation diagram of Eqs. (J3| in the (S,T) plane 
according to ([8]). A system with balance index e oscillates for 
E < S e (T) (hatched areas). Black dots indicate oscillating 
parameter sets of ©, with 6,8 e [10"\ 10], 6/a G [10, 1000], 
AG [1, 10 3 ], A/jLt G [10, 10 3 ], assuming the degradation mecha- 
nisms @ with a G [10 -6 , 1], k, x G [1, 100]. T c = 1 is the time 
scale at which transcriptional dynamics cannot be neglected, 
T op t = 2v2 is the location of the maximum of £i(T). 



a self-regulated gene with dynamical transcription rate: 
this circuit bifurcates most easily to periodic behavior, 
or more generally is least stable, at a finite value of the 
gene relaxation time given by r op t = 2\/2t c where 
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(10) 



The quantity r c gives the time scale at which dynami- 
cal behavior departs from the fast regulation case. Be- 
cause g<g is determined by ((3]), computing precise lower 
bounds on r c with (j 1 Of) requires specifying the degrada- 
tion mechanisms. Fixing Ao = (3q = 10 nm _1 , Po = 100, 
and g* = 0.5, provides an estimate t c = t c /Sm = 2 mn 
which is not unrealistically larger than typical gene in- 
duction times. At resonance, the oscillation period is 
T sc = 47r-\/2/3 r c (about 20 mn in the example above). 

Since a common interpretation of the diagram of Fig.[T] 
is that there are only two gene states (bound or unbound) 
(26| , one may wonder whether our deterministic analysis 
is relevant. If g is viewed as a temporal average of gene 
activity, our results are valid when the response time r is 
small compared to the oscillation period so that there are 
many binding/unbinding events by cycle 27] . Moreover, 
transcription is a complex process involving a number of 
distinct steps {4], and Eq. (fl"a|) is the simplest way to 
model memory effects arising from cooperativity in the 
transcription machinery Q. 

Anyhow, we now show that even when G is a stochastic 
variable jumping between and 1, our main result still 
holds: there is a time scale near T opt at which oscillations 
are enhanced. To this end we have carried out stochastic 
simulations of the reaction network of Fig. [1] using the 
Gillespie algorithm [28], varying the response time r at 
fixed Pq- Instead of regular oscillations, a sequence of 
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than one limiting step. 
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FIG. 3: Coefficient of variation of interpeak time intervals 
77 vs. mean residence time r for a stochastic simulation of 
© + ([9]) rescaled to original variables G, P, M at three values 
of threshold P (from top to bottom, P = 100, 500, 2500). 
Parameter values correspond to 5 = 1, = 1, a = d/Po, 
A = 21.54, n = 0.085, x = 95.5, K = 21.68, a = 10" 6 . n and 
T2 are the boundaries of the deterministic oscillation domain. 



irregularly spaced peaks in protein concentration is ob- 
served. A natural question is then whether protein peaks 
occur more regularly at parameter values where the de- 
terministic model oscillates, in particular when r = r opt . 

We define interpeak times as the time intervals At be- 
tween two crossings of P — 1.2P avg separated by at least 
one crossing of P — 0.8P avg , with P avg the mean pro- 
tein level (thus imposing a minimum amplitude of 40%). 
Their distribution is characterized by the coefficient of 
variation 77 = J 7 ^^ ■ A typical variation of rj with r 



<At>- 

in our system is shown in Fig. [3J It definitely suggests 
that the deterministic analysis remains relevant in the 
stochastic regime, since there is clearly a time scale near 
T opt where interpeak time fluctuations are minimal. 

In conclusion, we have shown that a nontrivial tran- 
scriptional dynamics can destabilize a self-regulated gene. 
Although it is known that nonlinear degradation mecha- 
nisms can induce oscillations in this system, we observe 
a resonance-like effect such that a much weaker nonlin- 
earity is required when the gene response time matches 
a characteristic time. Its expression can be computed 
analytically, which allows us to identify the parameter 
regions where this effect cannot be neglected. Stochas- 
tic simulations confirm the relevance of this time scale in 
the dynamics of the self-regulated gene. This shows that 
transcriptional dynamics is a possible source of oscilla- 
tor y b ehavior besides other deterministic [EBB 13 17, 
3, 13 13 El and stochastic [3 El effects. 

A natural question is whether our conclusions remain 
valid when more detailed transcriptional mechanisms or 
multiple sources of delay are taken into account. It can 
be shown that adding an explicit delay to Eqs. J3J only 
further destabilizes the circuit and thus does not essen- 
tially interfere with the effect of a gene response delay. 
Future studies should thus focus on cases where the tran- 
scriptional dynamics is more complex and features more 
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